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In the BCS theory of superconductivity, one assumes that all Cooper pairs have the same center of 
mass momentum. This is indeed enforced by self consistency, if the pairing interaction is momentum 
independent. Here, we show that for an attractive nearest neighbor interaction, this is different. 
In this case, stable solutions with pairs with momenta q and — q coexist and, for a sufficiently 
strong interaction, one of these states becomes the groundstate of the superconductor. This finite- 
momentum pairing state is accompanied by a charge order with wave vector 2q. For a weak pairing 
interaction, the groundstate is a d-wave superconductor. 
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In the original formulation of the BCS theory of super- 
conductivity [l| , all Cooper pairs are assumed to have the 
same center of mass momentum q. One possible gen- 
eralization of this theory is to introduce a pair ampli- 
tude for each center of mass momentum separately. In 
BCS theory for conventional superconductors, only one 
of these order parameters (OPs) is selected and the sta- 
ble state is the one where all pairs have the same mo- 
mentum and form the BCS condensate. For films in an 
external magnetic field, Fulde and Ferrell and, indepen- 
dently, Larkin and Ovchinnikov @, Q introduced a super- 
conducting (SC) state with coexisting pair momenta d=q, 
a state that explicitly breaks time inversion symmetry. 
For unconventional pairing symmetries, the competition 
between pair momenta is more complex and it has re- 
mained unresolved whether a bulk SC groundstate with 
different pair momenta may exist without magnetic field 

A SC state with different coexisting pair momenta gen- 
erally exhibits a spatially inhomogeneous charge density. 
One example of a superconductor of this type is the re- 
cently proposed "pair density wave" (PDW) state 0,0,0- 
It is characterized in real space by a two-component or- 
der parameter A(r) = A q e* qr +A_ q e _2qr . This structure 
bears some resemblance to the Larkin- Ovchinnikov state, 
but it preserves time inversion symmetry. The PDW is 
accompanied by a charge density pattern with wave vec- 
tor 2q. For this reason the PDW state has been proposed 
to describe the SC state of high-T c cuprates with coex- 
isting stripe order, especially Nd-doped La2-^Sr^Cu04 
and La 2 _ a 3a a; Cu04 for x = 1/8 jjj In 
particular, the recent experiments on the 1/8-doped ma- 
terial stimulated further theoretical studies to resolve the 
nature and the origin of the SC state in the charge or- 
dered phase [l3[. The PDW might be a candidate state, 
but so far a microscopic model which yields the PDW as 
its groundstate is lacking. 

In this article, we formulate an extended version of 
BCS theory using Gor'kov's equations for the Green's 



functions in the SC state. We explicitly allow for the 
coexistence of different finite-momentum pairing ampli- 
tudes in the absence of an external magnetic field. We 
identify conditions for a groundstate solution with finite 
OPs for the pair momenta ±q. This pairing state is re- 
alized beyond a critical interaction strength V c for an 
attractive nearest-neighbor interaction, and it is charac- 
terized by a charge stripe order, a gapless density of states 
(DOS) and a partially reconstructed Fermi surface. On 
the other hand, for V < V c , the d-wave superconductor 
is the stable groundstate. 

We start from a tight binding Hamiltonian on a square 
lattice with TV sites and periodic boundary conditions 



h J2 J2 V ^ k '' q)4 S C lk+q S ' C -k'+qs'Ck's. (1) 
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With a nearest and next-nearest neighbor hopping ampli- 
tude t and respectively, the single-electron dispersion 
has the form 



£k = — 2t [cos k x + cos k y ] + 4t' cos k x cos k y — ji 



(2) 



where \i is the chemical potential. 

For the superconducting state with singlet pair- 
ing, we use the BCS type mean-field decoupling 
scheme and approximate (c k ^ct k+ |C_k'+ q j,Ck't) — > 

(c k |C^ k+q |)(c_k / + q |Ck / t). The system is then repre- 
sented by the spin independent imaginary time Green's 
function C/(k, k',r) = — (T r c ks (r)c k/s (0)), and the 
anomalous propagators ^(k, k',r) = (T r c ks (Y)c_ k / s /(0)) 
and ^*(k,k',r) = (T r cl ks (r)cl sf (0)) for s ^ s f . The 
Heisenberg equations of motion for the normal and 
anomalous Green's functions lead to the Gor'kov equa- 
tions 1141: 
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£(k,k',cj n ) = g (k,u n ) 5 kk / - ^ A(k,q)f*(k - q,k',cj n ) 

q 

J^(k, k', cj n ) = <?o (k, c^ n ) ^ A(k, q)<?(-k', -k + q, -o; n ), 

q 

I 



(3) 
(4) 



where £7o(k, c<; n ) = [io; n — £k] _1 is the Green's function 
in the normal state and uo n = (2n — l)7rT is the fermion 
Matsubara frequency for temperature T. The order pa- 
rameter A(k, q) is determined by the self-consisteny con- 
dition 

A(k, q) = ~ E V ^ k '> k ' - <*' ( 5 ) 

n k' 

For the interaction, we choose a simple ansatz that al- 
lows for unconventional pairing; we assume an attrac- 
tive interaction between electrons on neighboring sites. 
In the singlet channel this is equivalent to the interac- 
tion V(k, k',q) = Vs(k, k',q) + Vd(k, k',q) in momen- 
tum space, with factorizable extended s-wave and d-wave 
components Vs(k, k',q) and V^(k, k',q), where 

V a , d (k, k', q) = Vg 8td (k - q/2)# M (k' - q/2). (6) 

V > is the attractive pairing interaction strength and 
<7s(k) = cos fc^ + cos fc y and #d(k) = cos fc^ — cos Thus 

A(k, q) = A s (q)# s (k - q/2) + A d (q)<fc(k - q/2). (7) 

The vector q labels mean-field solutions which corre- 
spond to order parameters in real space with phase wind- 
ing numbers q x and q y in x- and ^/-direction, respectively. 

If A(k, q) 7^ for a single momentum q / 0, then 
^-*(k, k',cj n ) and ^*(k, k', uo n ) have off-diagonal terms 
in momentum space, but (?(k, k',o; n ) is still diagonal. 
If A(k, q) 7^ for at least two different momenta q, 
then also <7(k, k', cj n ) has off-diagonal terms and the dis- 
crete translational invariance is broken. The charge den- 
sity is obtained from p{r) = l/iV^k k , e ir "( k_k/) n(k, k'), 
where n(k,k') = 2T £ n 0(k, k', w n ). Thus there 
are charge modulations whenever £/(k, k',a; n ) has off- 
diagonal terms. 

Inserting Eq. (|4|) into Eq. ([3]) leads to a system of cou- 
pled equations for the Green's function C/(k, k',a; n ). As- 
suming that n(k, k') <C n(k) for k / k', which will be 
verified a posteriori, ^(k, k + q, uo n ) and JF* (k, k + q, uo n ) 
are approximated by keeping only the term proportional 
to C/(k+q, k+q, uo n ) in Eq. (|4]) and by neglecting the sum 
over q otherwise. Within this approximation the Gor'kov 
equations are solved analytically, as shown below. 

In an ansatz for a self consistent solution of the 
Gor'kov equations (|3|) and (|4j), we choose Q trial vec- 
tors qi, . . . , qQ and set A(k, q) = for all other values 



of q 7^ q^. Thereby we test selected combinations of q- 
vectors for self consistent solutions. With this ansatz, the 
energy spectrum of the system consists of Q + 1 bands 
E a (k), a = 0, . . . , Q. The conventional BCS solution 
is realized for Q = 1 with just two quasiparticle bands 
and q = 0. Generally, one obtains a set of 2Q coupled 
self-consistency equations for A s (q^) and A^(q^): 



A M (q) 
V 



T 

~N 



^# M (k' - q^/2) ^V(k,k - q*,a; n ). 

k' n 

(8) 



Standard BCS theory provides isotropic solutions of 
Eq. (|HD with A s (0) = and A d (0) ^ for arbitrar- 
ily weak interaction strength V. Most remarkably, if 
V exceeds a certain interaction strength, we identify 
solutions with Q = 2 and with the specific set of q- 
vectors 0), (— q, 0)}. These solutions are anisotropic 
and have a subdominant extended s-wave contribution 
A s (q), which increases with increasing q. Below we will 
discuss in particular the time inversion symmetric zero- 
current solutions of Eq. (|5J), i.e., A s ^(q) = A s ^(— q). 

To test the stability of this solution, we solved Eq. (jHj) 
iteratively for selected combinations of q vectors and 
different initial values of the corresponding OPs. In 
particular, we investigated the stability of the above 
solution against decay into the q = state by us- 
ing the ansatz with the three center of mass momenta 
{((Z? 0)? ( — Qi 0) 7 (0? 0)}- We find that for a pure on-site 
interaction V(k, k 7 ,q) = Vb (s-wave pairing), finite mo- 
mentum pairing is unstable. All OPs with different q 
compete, even the ones with q and — q. Thus states with 
Q > 2 will always decay into a state with only one finite 
order parameter for on-site s-wave pairing. However, for 
the nearest-neighbor interaction (|6]), additional stable so- 
lutions emerge. The finite momentum pairing solutions 
are typically stable for a wide range of q values. Here 
the OPs for ±q do not compete, but rather support each 
other. The range of stability however decreases with de- 
creasing V, and eventually disappears. 

So far we have verified that stable finite momentum 
pairing solutions of the self-consistency equation exist. 
They refer to local minima of the free energy. To de- 
termine the groundstate at T = 0, the global minimum 
of the energy E = (H) = E k ^n(k) + EJA^q*) + 
A^(qi)]/V has to be determined with respect to all q. 
Figure [U shows the typical ^-dependence of E with a 
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FIG. 1: Energy E = (TL) as a function of pair momentum 
Q — (q, 0) f° r different pairing interaction strengths V. Cal- 
culations were performed for a 384 x 384 lattice with fixed 
electron density p — 0.8 and t' — 0.3t. For these parameters, 
the finite momentum pairing state becomes the groundstate 
for V > V c « 2.2t with q tt/3. 




FIG. 2: Density of states D(E) of the groundstate solutions 
for interaction strength V = 2t and V = 2.2t, corresponding 
to q = and q = respectively. The other parameters 

are the same as in Fig. [1] The coherence peaks of the q = 
state are split due to the van Hove singularity of the 2D tight- 
binding dispersion. 



minimum at q = and a further minimum for q > 0. 
The minimum at q = corresponds to the standard d- 
wave SC state. With increasing V, the energy of the 
minimum at finite q decreases accompanied by a shift to 
larger q. This implies the existence of a critical interac- 
tion strength V c , which depends on t' and the electron 
density p = 1/N n(k). Above V c the finite momen- 
tum pairing state is the groundstate. The optimal q sen- 
sitively depends on V, t' and p, but it is typically found in 
between q w tt/S and q « tt/2 for a wide parameter range. 
For even larger values of V, our calculations identified the 
groundstate to be of a checkerboard type with Q = 4 and 
with pair momenta {(g, 0), (— g, 0), (0, q), (0, —q)}. 

For the finite-g groundstate solutions the charge den- 
sity p(r) has an oscillatory part arising from the off- 
diagonal terms of the Green's function. For the q = 
(=bg, 0) state the charge density forms a sinusoidal stripe 
pattern with wave number 2q. Correspondingly, the 
charge density varies as 



p(r) = p + pi cos(2gx), 



(9) 



with an amplitude pi/ p « 2%, which justifies the as- 
sumption of small charge modulations in the above ap- 
proximation for ^-*(k, k',o; n ). For q = 7r/3, the wave- 
length of the stripe pattern is therefore three lattice con- 
stants. The charge modulation in the SC state suggests 
to include a self-consistent charge density- wave (CDW) 
OP in the mean-field decoupling scheme of the Hamilto- 
nian (pQ). We have analyzed this extension with coexist- 
ing OPs for SC and CDW order for selected cases. The 
CDW OP tends to stabilize the state with finite momen- 
tum pairing but it remains small and does not change 
the solutions qualitatively. 

The finite momentum pairing state has further char- 
acteristic properties that are at variance with a BCS 



like d-wave superconductor (with q = 0). The DOS 
D(E) = X]kl m( ^(k>k> E — i0+), where ImGis the imag- 
inary part of the analytical continuation of Q to the real 
frequency axis, is shown in Fig. O For q = 7r/3, the 
DOS bears little resemblance to a d-wave like gap as the 
coherence peaks are split and the DOS is finite at the 
Fermi energy. A similar splitting is observed for current 
carrying d-wave states 15|, [l6| which originates from the 
Doppler shift of the finite momentum eigenstates. 

Figure [3] displays the characteristic momentum space 
properties of the finite momentum pairing state and, for 
comparison, of the q = d-wave superconductor. In the 
finite momentum pairing state, the momentum distribu- 
tion function n(k) = n(k, k) develops structures with 
sharp boundaries. These boundaries consist of lines in 
momentum space with E a (k) = 0, for a = 0, 1 or 2 in- 
dicative of a Fermi surface reconstruction. The zero en- 
ergy states generate Fermi- arc like structures as shown in 
Fig- Eli. For t' = 0, n(k) is similar to the result obtained 
in Ref. @. The pair density P(k) = £ • P(k, k - q*), 
where P 2 (k, k 7 ) = 2(cj c |cl k ,|)(c_j c ^Ckt), clearly shows 
the fingershaped k-space structures of n(k) contain un- 
paired electrons only. The overall number of pairs is 
smaller in the finite momentum pairing state than in the 
q = state. This seems to contradict the fact that it has 
the lower energy. The latter, however, consists of both 
the kinetic energy, which rises in the SC state and acts 
against the formation of pairs, and the gain of conden- 
sation energy. The optimal balance between these two 
contributions depends on details of the single particle ki- 
netic energy and the interaction potential V(k, k 7 ,q) 
and does not generally favor a larger number of paired 
electrons. 

In this article we have shown that the extended BCS 
theory with attractive nearest neighbor interaction pro- 
vides self consistent solutions with the simultaneous for- 




FIG. 3: Momentum space properties of the finite momen- 
tum pairing state with q = tv/3 and the same parameters as 
in Fig. [2] (right panels) and for comparison the d-wave su- 
perconductor for q = (left panels), (a), (b): Occupation 
probability function n(k). (c), (d): Density of states with 
zero energy ImC(k, k, — id) (here: S = 0.04£). (e), (f): Pair 
density -P(k). 

mation of electron pairs with center of mass momenta q 
and — q. It is a microscopic solution which constitutes 
a stable macroscopic state of the PDW type which was 
proposed to describe the striped SC phase in hole doped 
214 cuprates. This finite momentum pairing state is the 
groundstate beyond a critical interaction strength V c . V c 
depends sensitively on the band filling p and ranges from 
V c « I At for p = 0.6 to V c « 3M for p = 1. This is con- 
sistent with the result in Ref. Q that only the uniform 



phase with fixed q can be the groundstate of the BCS 
Hamiltonian in the weak coupling limit. 

The results described above as solutions of Gor'kov's 
equations are alternatively obtained by diagonalizing the 
mean- filed decoupled Hamiltonian (pQ) using a Bogoliubov 
transformation or by solving an extended version of Bo- 
goliubov - de Gennes equations in real space [l7| . 

Our results demonstrate as a proof of principle that 
stable groundstate solutions of the pairing Hamiltonian 
(PP) exist with coexisting finite momentum pairing am- 
plitudes for center of mass momenta q = 0) and 
— q; these solutions are absent for an attractive con- 
tact interaction. Due to the concomitant striped charge 
density modulation with wave vector 2q, a connection 
to the striped superconductor La 15 / 8 Ba 1 / 8 Cu04 appears 
tempting. However, without the inclusion of additional 
correlation effects as the source for a possible spin order 
pattern, we consider it premature to draw conclusions 
about the favorable wavelength of the stripes. 
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